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Abstract 



O 

. We investigate the thermodynamic behaviour of a Bose gas interacting 

\^ 

On ■ with repulsive forces and confined in a harmonic anisotropic trap. We de- 



velop the formalism of mean field theory for non uniform systems at finite 



. temperature, based on the generalization of Bogoliubov theory for uniform 

O ■ gases. By employing the WKB semiclassical approximation for the excited 



states we derive systematic results for the temperature dependence of vari- 
ous thermodynamic quantities: condensate fraction, density profiles, thermal 
energy, specific heat and moment of inertia. Our analysis points out impor- 
tant differences with respect to the thermodynamic behaviour of uniform Bose 
gases. This is mainly the consequence of a major role played by single particle 
states at the boundary of the condensate. We find that the thermal depletion 
of the condensate is strongly enhanced by the presence of repulsive interactions 
and that the critical temperature is decreased with respect to the predictions 
of the non- interacting model. Our work points out an important scaling be- 
haviour exhibited by the system in large N limit. Scaling permits to express all 
the relevant thermodynamic quantities in terms of only two parameters: the 



reduced temperature t = T/T^ and the ratio between the T = value of the 
chemical potential and the critical temperature for Bose-Einstein conden- 
sation. Comparisons with first experimental results and ab-initio calculations 

are presented. 
02.70.Lq, 67.40.Db 
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I. INTRODUCTION 



The experimental realization of Bose-Einstein condensation (BEC) in atomic gases |l| 
confined in magnetic traps has stimulated a novel interest in the theory of interacting Bose 
gases. In particular the recent experiments carried out at Mit on ^^Na and at Jila on 
^^Rb have provided first quantitative results for relevant thermodynamic quantities, such 
as the temperature dependence of the condensate fraction and the release energy. Based on 
the pioneering works by Bogoliubov and Huang, Yang, Luttinger and Lee the theory 
of interacting Bose gases has so far been mainly focused on uniform systems. The possibility 
of extending it to non-uniform gases, first explored in the works by Gross, Pitaevskii and 
Fetter and stimulated by the recent experimental activity, is expected to open new 

challenging perspectives in many-body physics, since these systems are characterized by 
new length scales and symmetries. The thermodynamic behaviour of trapped Bose gases 
has been already the object of several theoretical works, starting from the investigation of 
the ideal Bose gas §], the development of the Hartree-Fock formalism the inclusion 



of interaction effects using the local density approximation [|TO|,|TT|, up to the most recent 



approaches based on self-consistent mean-field theory [T^-|15| and numerical simulations [jT6 

The purpose of this paper is to provide a systematic discussion of the properties of these 
systems at thermodynamic equilibrium. We will always assume that the dilute gas condition, 

na^ « 1 , (1) 

is satisfied. This condition is reached in all the available experiments. Here n is the atomic 
density whose typical values range from 10^^ to 10^^ cm~'^, while the triplet spin s-wave 
scattering length a can reach a few tens Angstrom. The corresponding value of na^ is 
then always extremely small (10~^ - 10~^) and one would naively expect that the role of 
interactions is a minor effect in such systems, at least for the equilibrium properties. This is 
not the case for these trapped atoms, because of the combined effect of trapping and Bose- 
Einstein condensation. To illustrate this important feature let us consider the simplest case 
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of N atoms occupying the ground state of an harmonic oscillator trap. The ratio between 
the interaction energy and the quantum oscillator energy hu behaves as: 

~ , (2) 

JMnuj ano 

where ano = {h/muy^'^ is the harmonic oscillator length. In the available experiments 
O'/O'HO is ~ 10^'^ while ranges between 10^ and 10^ so that the ratio @) can be very 
large, revealing that at least the ground state properties can never be calculated treating the 
interaction potential in a perturbative way. Differently from what happens in a homogeneous 
Bose gas where the density is kept fixed, the repulsion has here the effect of expanding the 
condensate wave function. Actually, for sufficiently large values of N , it happens that 
the T = equilibrium properties of these new systems are governed by the competition 
between the 2-body force and the external potential, the kinetic energy playing a minor 
role. In this regime the ground state density is well approximated by the Thomas-Fermi 
(TF) approximation 

Tfl 

^^^^ = rT2-(/^0 - Ve^t{r)) , (3) 

if VJ=Kt(i") < A^o and n(r) = elsewhere. Here /io is the T = chemical potential, fixed by the 
normalization condition. Assuming a harmonic potential V^xt = Tn{ijj1x'^ + ^yU"^ + t^2-2^)/2, 
the chemical potential takes the useful expression [p!7tp!8| 



/zo = — ISAT , (4) 

z \ aHoy 

where the frequency uj is the geometrical average 

u = [uxUJyUj^y^^ (5) 
of the oscillator frequencies and fixes the harmonic oscillator length 

The root mean square (r.m.s.) radius of the condensate takes the form 



= (^15iV—j (7) 

where A is the anisotropy parameter of a cyhndrically simmetric trap A = io^j^z = ^y/^z- 
In the presently available traps one has typically ano = 10~^ cm and J (r^) ranges from 2 



to 100 xlO~^ cm. The validity of the TF approximation (|^), (|) is based on the condition 
Na/auo ^ 1 which also implies fio/hu ^ 1. These results have been obtained at T = 0. 
They play however a relevant role also at finite temperature. In particular the value of 
/io fixes an important reference in the scale of temperatures. For fc^T smaller than /xq the 
thermodynamic behaviour of the system is in fact affected by collective effects (phonon 
regime in a homogeneous superfluid). Conversely for higher temperatures it is dominated 
by single-particle effects. 

It is highly interesting to compare the value of the chemical potential /iq with the one 
of the critical temperature. Due to the minor role played by interactions at high T, a 
reliable estimate of the critical temperature is obtained using the non-interacting model. In 
this model the Hamiltonian is given by if = Y^i{pf/2m + Vexti^i)) and its thermodynamic 
behaviour can be worked out exactly. For sufficiently large numbers of atoms this model 



predicts the value [|8|,P0 



ksT^ = j (8) 

for the critical temperature. In the presently available traps one has typically huj/ks = 3 — 5 
nK, while T° ranges from 50 to 300 nK. Using the TF result (^) the ratio r] between the 
chemical potential and the critical temperature takes the useful form 

= = a{N'/'^f"> ^ 2.2A{a'nT=o{r = 0))'/' . (9) 

ksT^ ano 

where a = |C(3)^/^15^/^ ^ 1.57 and we have used result (0) for the density. This equation 
reveals an interesting behaviour of these trapped Bose gases. In fact although the gas 
parameter na^ is always very small, its 1/6-th power can be easily of the order of unity 
and in actual experiments it turns out that /io is 0.3-0.4 /c_bT°. Furthermore the ratio 
depends on the number of atoms in the trap in an extremely smooth way (rj ~ N^^^^). 



The ratio ^Q/ksT^ plays a crucial role to characterize the interplay between interaction 
and thermal effects. In particular we will show that the thermodynamic properties of the 
Bose condensed gases, for sufficiently large values of A^, depend on the relevant parameters of 
the system (number of atoms in the trap, s-wave scattering length a, oscillator length and 
deformation of the trap) through the dimensionless scaling parameters rj. Note that such a 
parameter involves a different combination of the quantities N and ajauo-, as compared to 
the ratio @ which is instead the natural scaling parameter for the ground state properties 

Mil- 



In the paper we will investigate the thermodynamic behaviour of the trapped Bose gas 
using a mean-field theory based on already available extensions of Bogoliubov theory to 
finite temperature. In particular we will work in the so called Popov approximation. The 
details of this approximation are discussed in Sect. 2. Here we recall some major features. 
At T = the theory exactly accounts for the Bogoliubov spectrum of elementary excita- 
tions. At high T it approaches the finite-temperature Hartree-Fock theory with which it 
exactly coincides above the critical temperature. These two important features (collective 
and single-particle effects) are the main ingredients of the Popov approximation, which is 
then expected to be a good approximation in the whole range of temperature except near 
where mean-field theories are known to fail. An interesting feature of these trapped gases, 
which distinguishes them from uniform Bose gases, is that single-particle excitations play an 
important role also at low temperatures. The reason is that these systems exhibit two kinds 
of low lying excitations: on the one hand one has collective modes of the condensate, on 
the other hand one has single-particle excitations near the classical boundary. In the paper 
we discuss the relative importance of the corresponding contributions to thermodynamics. 
We find that the single-particle contribution is always relevant. In particular, even at low 
temperatures, the thermal energy is governed by the excitation of single-particle states dif- 
ferently from what happens in homogeneous Bose systems. This fact explains why in the 
trapped gases the Hartree-Fock approximation provides a very good description for most 
thermodynamic quantities in a wide range of temperatures. 



In order to simplify the formalism and to carry out explicit calculations in a systematic 
way we will restrict our analysis to temperatures such that 

ksT » hu . (10) 

This assumption allows us to describe the motion of the thermally excited atoms in the WKB 
semi-classical approximation In a deformed trap condition ( p!o|) is separately imposed 



on the three oscillator energies hux, huy and hujz- Current experimental investigations 
are carried out at temperatures well consistent with the above condition. The study of the 
thermodynamic behaviour for lower temperatures (of the order of the oscillator temperature) 
requires a quantum description of the low lying elementary modes, beyond the semi-classical 
approximation. These modes have been the object of recent extensive experimental and 



theoretical work [22-25.14 



The final equations resulting from the mean-field approach have the form of coupled 
equations for the order parameter and for the density of the thermally excited atoms. Their 
numerical solution will be systematically discussed during the paper and compared with 
the available experimental data as well as with the first results obtained with Path Integral 
Monte Carlo simulations with which we find a remarkable agreement ||T^|. 



The paper is organized as follows. Sect. 2 is devoted to the presentation of the theoretical 
formalism. We discuss the Popov approximation, we obtain formal results for the relevant 
thermodynamic quantities and we derive explicitly the scaling behaviour as well as the low 
T expansion of the thermodynamic functions. In Sect. 3 we discuss our results for the 
temperature dependence of the condensate fraction, density profile, energy, specific heat 
and moment of inertia. Finally in Sect. 4 we make a comparison between the properties of 
a trapped Bose gas and the ones of the uniform gas. 



II. THEORY 
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A. Self-consistent Popov approximation 



The self-consistent mean-field method described here follows the treatment of the inho- 
mogeneous dilute Bose gas developed in Ref. [0] within the Bogoliubov approximation at 
T = and its extension to the finite-temperature Hartree-Fock-Bogoliubov (HFB) approxi- 
mation recently discussed by Griffin in Ref. |r^|. The starting point is the grand-canonical 



Hamiltonian written in terms of the creation and annihilation particle field operators ^/'^(r) 
and '?/'(r) 

+ I y rfrV^t(r)7/>t(r)7/^(r)V^(r) . (11) 

In the above equation Vext{T^) = m{ujlx'^+u!yy'^+uj'^z^) /2 is the anisotropic harmonic potential 
which provides the confinement of the atoms and g = Anfi^a/m is the interaction coupling 
constant fixed by the s-wave scattering length a. By separating out the condensate part of 
the particle field operator one can write ipi^) fhe sum of the spatially varying condensate 
wave-function $(r) and the operator '?/'(r), which acts on the non-condensate particles, 

^(r) = $(r) + V^(r) . (12) 

The wave-function $ is formally defined as the statistical average of the particle field op- 
erator, $(r) = (■?/'(r)), and represents the order parameter of the system in the condensate 
phase where gauge symmetry is broken |2^. Using decomposition (|1^) the interaction term 
in (p!T| ) becomes (all quantities depend on r) 

^t^t^^ = |$|4 + 2|$p$> + 2\^\'^^ip^ 

+ Al^'^ip^ + ^*^'ipip + ^^ip^ip^ (13) 

In the following the cubic and quartic terms in ■?/', ip'^ in the last line of the above equation 
will be treated in the mean-field approximation. For the cubic terms one gets 
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'4;\r)tp\r)tP{r) = 2{4j\r)4j{r))'^\r) = 2nT{r)4j\r) , (14) 

where nT{r) is the non-condensate density of particles. Whereas the quartic term becomes 

^t(j.)^"^(r)^(r)^(r) = 4(^'^(r)^(r))^"^(r)V^(r) = 4nT(r)^'^(r)^(r) . (15) 

To obtain the mean-field factorizations (|T^) and (|I^) we have neglected the terms propor- 
tional to the anomalous non-condensate density mT{r) = {tp{r)ilj{r)), and to its complex 
conjugate. This approximation, which is discussed in depth in Ref. [Q, corresponds to the 



so called Popov approximation and is expected to be appropriate both at high temper- 
atures, where ht ^ ttlt, and at low temperatures where the two densities are of the same 
order, but both negligibly small for the very dilute systems we are considering here. Us- 
ing decomposition (|12D and the mean-field prescriptions ([T^), and ([TSD the grand-canonical 



Hamiltonian (|TI]) can be written as the sum of four terms 

K = Ko + Ki + Kl + K2 . (16) 

The first term Kq contains only the condensate wave function $(r) 

Ko = J rfr$*(r) (^-^ + V,^t{r) - /i + |no(r)^ <l>(r) , (17) 

where we have introduced the condensate density no(r) = |$(r)p. The Ki and kI terms 
are linear in the operators ip, ip'^ 

K, = J dvij\v) [-^ + Ve^tij) + g{no{r) + 2nT{r))j $(r) , (18) 

and finally K2 is quadratic in ip, ip^ 

K2 = J dr'4;\r)Cip{r) + dr^^ {r)ip\r)ip\r) 



+ dr^*\r)^{r)^{r) . (19) 
In eq. (|T^, to make notations simpler, we have introduced the hermitian operator 



C = -— + Kxi(r) + 2gn{v) 
2m 



(20) 



where 



n(r) = no(r) + nxir) 



(21) 



is the total particle density of the system. The Ki, K\ terms in the sum (|T6]) vanish if $(r) 
is the solution of the equation for the condensate 



-7^ \- Vextij) - + 5'(^o(r) + 2nT(r)) 



<l>(r) = [C- gn^iv)] <l>(r) = , 



(22) 



with C defined in eq. (|20|), whereas the quadratic term K2 can be diagonalized by means of 
the Bogoliubov linear transformation for the particle operators ip, ip'^ 



(23) 



In the above transformation a^, and a] are quasi-particle annihilation and creation operators 
which satisfy the usual Bose commutation relations, whereas the functions %(r), fj(r) obey 
to the following normalization condition 



(24) 



The grand-canonical Hamiltonian (|Tl|) , in the Popov mean-field approximation, takes the 
form 



Kpopov = J dr^r) + Ve.t{r) - /i + |^o(r)^ $(r) + ^ e^a] 



a 



(25) 



where the real wave-function $(r) is the solution of eq. ( p2D and the quasi-particle energies 
tj are obtained from the solutions of the coupled Bogoliubov equations 



Cvj{r) + gno{r)uj{r) = -ejVj{r) 



(26) 
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for the functions Uj{r) and Vj{r). In eq. ( pS]) we have neglected the term — J2j I dr\vj{r)\'^ 
which at T = accounts for the energy of the non-condensate particles in the ground 
state. This is a good approximation since, for dilute systems satisfying condition (|T]), the 
quantum depletion J2j I dT\vj{r)\^ of the condensate is a minor effect. Explicit calculations 
show that this depletion is always smaller than one percent of the total number of atoms 
in the condensate IQ. We have therefore at T = 0, ip{r) = $(r) and nT{r) = mT{r) = 0, 



whereas at finite temperatures we take ht and proportional to the thermal density of 
quasi-particles 



riT 



j:i\u,\' + \v,\'){a]a,) , (27) 



and 



rriT = 2^^UjV* (a'jaj) . (28) 
j 

Equations (|26| ) have the typical form of the equations of the random phase approximation 



and have been recently solved in Ref. at T = 0; the resulting energies e,- give the 



excitation frequecies of the collective modes of the system. Very recently, first solutions 



have become available also at finite temperature WM. The coupled equations (EH) can be 



solved very easily if one uses the semiclassical ||13[ or local density approximation |110| , p!l 
We write 

Mj(r) = m(p, r)e''^W Vj{r) = v{p, r)e^^('-) , (29) 

where the impulse of the elementary excitation is fixed by the gradient of the phase p = hVip 
and satisfies the condition p ^ Ti/ano- We also assume that ^(p, r), w(p, r) (and hence tlt) 
are smooth functions of r on length scales of the order of the harmonic oscillator length 
o-HO = {h/mujy^'^. Since the typical size of the thermally excited cloud is of the order of 
the classical thermal radius {2kBT/muj^y^'^, the semiclassical approximation makes sense 
only if ksT ^ hu |^T|. The sum over the quantum states j labelling the solutions of the 
Bogoliubov equations (^6|) is replaced by the integral / dp / {2nh)^ . By neglecting derivatives 
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of u and v and second derivatives of ip it is possible to rewrite the Bogoliubov-type equations 
(E6|) in the semiclassical form 



1^ + Vextir) - /i + 2gn{r) \ u{p, r) + gno{v)v{p, r) = e(p, r)u(p, r) 



^ + Vext{r) -11 + 2gn{Y) \ t;(p, r) + 5(no(r)M(p, r) = -e(p, r)w(p, r) 



(30) 



The solution of equations (|30| ) is easily obtained and one finds the following result for the 
spectrum of the elementary excitations 

e{p,v)={^^ + VUv)-^i + 2gn(v)^ - /^oW j , (31) 

with the coefficients m(p, r) and f (p, r) given by 

2 , , pV2"^ + Vext{Tc) - /i + 2gn{Y) + e(p, r) 
u (P,r) = 



2e(p,r) 

2( ^ _ P^l2m + Vextj'c) - fi + 2gn{r) - e(p, r) 

2e(p,r) 

«(p,r).(p,r) = -^-^ . 

The spectrum (|3T|) for the elementary excitations is significant only for energies e(p, r) ^ /icu. 
In fact the low-lying modes of the system are not properly described in the semiclassical 
approximation and their study requires the solution of the coupled Bogoliubov equations 



, P^ . As already mentioned above, the condition e(p, r) huj restricts the study 
of thermodynamics to temperatures much larger than the oscillator temperature huj/ks- In 
practice the resulting analysis is useful only for systems with large values of A^, where the 
critical temperature is much larger than the oscillator frequency (see eq.(|D). 

It is interesting to analyze in more detail the excitation energies (^l]). First of all in the 
absence of interparticle interactions {g = 0) the elementary excitations reduce to the usual 
single particle spectrum in the presence of the external potential Vext 

2 

e°(p,r) = |^ + V;,,(r)-/i , (33) 
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where the chemical potential is given by = 3^ci}/2, and uj = + Uy + u!z)/3 is the 
arithmetic average of the oscillator frequencies. In the limit Na/auo ^ 1 (Thomas- Fermi 
limit) we can distinguish two regions in space depending on whether the coordinate r lies 
inside or outside the condensate cloud. Inside the condensate region one can neglect the 
kinetic term in eq. (^) and one has approximately gijiQ + ^tit) = — V^xt- In this region 
eq. (|3T|) becomes 

^^ns^de^^^ r) = £- (p^ + Amgn,{v)) , (34) 
which corresponds to the local-density form of the Bogoliubov spectrum. For small p eq. 



(0) gives the phonon-like dispersion g«"««f^e(^p^ r) = c(r)p where c(r) = ygno{r)/m is the 
local velocity of sound. Far from the condensate cloud one instead has Vext » '^9 {no + tit) 
and we are left with 

6°"*^^'^(p,r) = |^ + \/e.t(r)-/x , (35) 

which has the same form as the independent particle spectrum (0), but now the chemical 
potential contains the effects of interaction and is in general much larger than its non- 
interacting value. Finally it is worth noting that when the energy e(p, r) is much larger 
than the chemical potential fi one obtains the Hartree-Fock type spectrum ||^ 

2 

e^^(p, r) = ^ + V,xt{r) - /i + 2gn{r) . (36) 
2m 



In the semiclassical scheme the canonical transformation (p3|) becomes the familiar Bo- 
goliubov transformation in momentum space relating the particle and quasi-particle creation 
and annihilation operators 

^(p, r) = m(p, r)a(p, r) + v{p, r)a\-p, r) , 

V^"^(p,r) = M(p,r)«1'(p,r) +t;(p,r)a(-p,r) . (37) 

The quasi-particle distribution function /(p, r) is obtained from the statistical average 
(a"''(p, r)a(p, r)), given by the usual Bose distribution 
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/(p, r) = (at(p, r)a(p, r)) = ^.^(^(p^ _ , , (38) 

whereas the particle distribution function -F(p, r) is given by the statistical average 
(^/^^(p, r)ip{p, r)) of the particle field operators in the coordinate-momentum representation. 
By making use of the transformation (R^) one gets 



F(p, r) = {ip^p, r)?/5(p, r)) = (^^(p, r) + v^{p, r) j /(p, r) 

/ de{p,r)\ 



/(P,r) , (39) 

no,nT 



with 



'9e(p,r)\ /2m + Vext{r) - jJ, + 2gn{r) 



da I e(p, r) 

" / no,nT ' 



(40) 



It is worth noticing that the particle and quasi-particle distribution functions differ only in 
the low-energy regime. In fact at high energies the excitation spectrum takes the HF form 
(^), de/dfi = —1 and hence F(p, r) = /(p,r). Conversely in the phonon regime one has 
F(p,r) = mc(r)/(p,r)/p. 

The non-condensate density nrp can be obtained in a self consistent way by integrating 
over momenta the particle distribution (p9|) 



and finally the chemical potential /i is fixed by the normalization condition 

N = No{T) + Nt = J rfmo(r) + J drnrir) , (42) 

where No{T) and Nt are respectively the total numbers of atoms in the condensate and out 
of the condensate in equilibrium at temperature T. 



Equations (0), ( |3TD and (|38|)-(^2D must be solved self-consistently. The solution of these 
equations yields directly, for a given temperature T, the chemical potential /i, the condensate 
density no{r) and the density of thermally excited atoms ^^(i')- 

Starting from the results discussed above it is possible to study the thermodynamic 
properties of the system whose excited states can be classified in terms of the states of a gas 
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of undamped quasi-particles with energy spectrum (pll). The entropy of the system is then 
readily obtained through the combinatorial formula ||28| 

One must notice that the dependence of the entropy on T is determined also through the 
temperature dependences of the quantities fi, Uq and which enter the excitation spectrum 
e(p, r) (see eq. (|3ll) ). Starting from the entropy one can work out all the thermodynamic 
quantities. For example the total energy of a A^-particle system at temperature T can be 
calculated using the thermodynamic relation 

E{T) = E{T = 0) + r dT' Cn{T') , (44) 

JO 

where 

C«m^T{f)^ (45) 

is the specific heat at constant N and E{T = 0) is the energy at zero temperature given by 
the well known Gross-Pitaevskii energy functional [^,0. 



E{T = 0) = I dr^r) (^-^ + 14.*(r) + |no(r)^ $(r) , 



(46) 



where $(r) is here the condensate wave function at T = 0. 

In the last part of this section we discuss the rotational properties of a trapped Bose gas 
which point out the effects of superfluidity p9|. The moment of inertia B, relative to the 



z direction, is defined as the linear response of the system to the perturbation field —QL^ 
according to the relation 

(L.) = no , (47) 

where is the 2;-component of the angular momentum. In terms of the particle field 
operator one has 

L, = J dr^^r) 4 ^(r) , (48) 
15 



where £2 = —ih{xd/dy — yd/dx). The statistical average in eq. (|^) must be calculated in 
the ensemble relative to the perturbed grand-canonical Hamiltonian K'{VL) = K ~ QL^. In 
the frame rotating with angular velocity Q the quasi-particles have energies 

en(p,r)= |^(^|- + V;,i(r)-/i + 2(7n(r)^ -/^oWj - f](r x p), , (49) 
and are distributed according to the Bose function 

/n(p,r) = . / ^ . (50) 

exp(ef7(p,r)) - 1 

In the limit of an axially symmetric trap [ux — »■ ujy) , the condensate wave function does not 
contribute to the average (L^) which is determined only by the angular momentum carried 
by the elementary excitations: 

= I (^^"""P^^^^^P'"^ • ^^^^ 
To lowest order in Q one obtains the following relevant result for the moment of inertia 

Eq. (|52D can be also written in the useful form 

6 = 1 dr{x^ + y^)pn{r) , (53) 

where 

provides a generalization of the most famous Landau formula for the density of the normal 
component to non homogeneous systems. The deviation of the moment of inertia from the 
rigid value 

Qrig ~ ^ J dr{x'^ + y'^)n{v) , (55) 

provides a "measure" of the superfiuid behaviour of the system (see section 3-C). 
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B. Hartree-Fock approximation 



In the previous section we have seen that outside the condensate and in general for high 
excitation energies the Popov quasi-particle spectrum (^) reduces to the famihar Hartree- 
Fock spectrum (|36|). The self-consistent HF method has been already employed to study 
the thermodynamic properties of inhomogeneous dilute Bose systems 0,0. 

The HF approximation can be easily obtained by neglecting the small " hole" components 
Vj in the first of the Bogoliubov equations (|2BD. The resulting equation then takes the form 



Cuj{r) = (^-^ + Vext{r) - /i + 2gn{r)j Uj{r) = ejUj{r) . (56) 

Notice that the Hamiltonian (|56|) differs from the effective Hamiltonian (p2D yielding the 
condensate wave function (for a more detailed derivation and discussion of the Hartree-Fock 
approximation for Bose gases see Ref. 0). By applying the semiclassical prodedure to solve 
eq. (|56[) , one finds the energy spectrum enrip,''^) of eq. (|36D. In the HF approximation 



there is no difference between particles and quasi-particles and therefore the distribution 
functions /(p, r) and F(p, r) coincide {denp/dfi = —1) 

F(p, r) = /(p, r) = , \ ^. ^ . (57) 

exp(e/^ir(p, r)/fcBJ ) - 1 

Results (P3|)-(^5D for the thermodynamic quantities obtained in section 2-A hold also in 
the HF approximation if one simply replaces the excitation spectrum e(p, r) with the cor- 
responding HF spectrum (|36D . By calculating the statistical average of the grand-canonical 
Hamiltonian ( pT]) over the eigenstates eq. (^) one obtains the following explicit expression 



for the energy of the system in the HF approximation 

+ J drVe.t{r)n{r) + ^ J dr{2n\r) - nl{r)) . (58) 

In the above equation the first two terms give respectively the kinetic contribution from 
the condensate and non-condensate particles. The third term gives the confining energy 
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averaged over the total density n, and the last term represents the interaction contribution. 
The HF energy functional (|58D is consistent with the thermodynamic relations ((^3|)-(^5|)). 

Results (|5^)-(|5^ for the moment of inertia also hold in the HF approximation with the 
distribution /(p, r) given by (0). In the HF approximation one easily finds p„(r) = mnT{r) 
and hence the moment of inertia is given by 

Qj^p = m j dr {x^ + y^) nr(r) (59) 

consistently with the finding of Ref. The identification of the normal component p„ 

with the non-condensate density is a peculiar property of the HF approximation. In the 
Popov approach these two quantities behave differently at low temperatures because of the 
presence of phonon-type excitations. 



C. The non- interacting model 

The results for the thermodynamic functions of the trapped ideal gas are easily obtained 
by setting the coupling constant g equal to zero in the equations derived above. The equation 
for the condensate (|2^ ) takes the form 

(-^ + K.*(r)-/ij<I'(r) = , (60) 

whose solution is the gaussian function 

$(r) = Y^iVo [^-^j exp i^- — {uj,x^ + ujyy^ + uj,z^)j , (61) 

normalized to the number of condensate atoms / (ir|$p = A^o, whereas the chemical potential 
is /i = h{u!x + ^^y + t^z)/2. The semiclassical spectrum ( pT]) takes the form (|33|) and the non- 
condensate density is given by 

nT(r) = / ^^(exp((pV2m + K.*(r)-/x)ABr)-l)-i 

= A^3^3/2(exp((V;xt(r) - fx)/kBT)) , (62) 
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where At = h^2TX /mkBT is the thermal wavelength and ga{z) = Z^^i -2^"/''^" is the usual 
Bose function. The total number of non-condensate atoms is obtained by integrating ut 
over space 



Nt = J dr nrir) = (^^^^ gsiexpifi/ksT)) . 



(63) 



In the region of temperatures where the semiclassical approximation for the excited states 
is valid, ksT ^ Tioj, one can expand the gz{,z) function around z = 1 and obtain to lowest 
order in N^^^^ 

N ^ ^2{a3)y/'uj ^ ' ^^^^ 

where we have introduced the reduced temperature t = T/T^, and T° = (N / ((3)y^^huj / ks 
is the critical temperature in the large limit. Finite size effects, which are accounted for 
by the second term in eq. (|64D , are proportional to the ratio u/uj between the arithmetic 
LJ = + Uy + ujz)/^ and the geometric averages of the oscillator frequencies and thus 
depend on the trap anisotropy parameter A = ujz/^x = ^z/^y, having a minimum for an 
isotropic trap (A = 1). These finite size effects shift the value of the critical temperature 



from T° towards lower temperatures. By setting the right hand side of eq. (|6^) equal to 



unity one finds the following result for the shift in the critical temperature |30 



TO 2(C(3))2/3u; uo ^ ' 

The temperature dependence of the energy per particle is easily obtained in the semiclassical 
approximation and one has (for t < 1) 

It is worth noticing that the shift ( |65D in the critical temperature as well as the A^~^/^ terms 
in eqs. (p^), ( |66D exactly coincide with the results obtained to lowest order in A^^^^'^ from a 



full quantum treatment of the non- interacting model in the large limit |^ . For = 1000 



the predictions (|64|), (pq) are already indistinguishable from the exact result obtained by 
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summing explicitly over the excited states of the harmonic oscillator Hamiltonian. This is 
well illustrated in Fig. 1 where we plot the prediction Nq/N (see eq. (|6^) for the tempera- 
ture dependence of the condensate fraction compared to the exact quantum calculation and 
the large N behaviour Nq/N = 1 — obtained by setting = in eq. (p3|). The accuracy 
of the semiclassical approximation reveals that the main origin of finite size effects in the 
thermodynamic properties of these systems arises from the quantum effects in the equation 
for the condensate, which are responsible for the proper value of the chemical potential, 
rather than from the quantum effects in the equation for the excited states. 

It is also important to notice that in the large N limit all the thermodynamic properties in 
the non-interacting model depend only on the reduced temperature t = T/T^ (see eqs. (|64D, 
(|66|) with oo). In the next section it will be shown that in the presence of interaction 
the thermodynamic properties of the system depend also on the additional parameter r) 
defined in eq. (|^). 

Finally, the ratio Q/Qrig of the moment of inertia to its rigid value has been calculated 
for the non-interacting model in Ref. ES]. In terms of the reduced temperature t this result 



reads 



1-^ ..s , -..4M.m/.s ' (67) 



Qr^g 1 - + jt^XNy/^ ' 

where 7 = 2C(4)/(C(3))^/3 ^ 1.69. Starting from (^) one finds that in the non-interacting 
model the effects of finite size are much more important for O than for the other thermody- 
namic quantities. 



D. Thomas- Fermi limit and scaling 

When the number of atoms in the condensate is enough large the Thomas-Fermi ap- 
proximation, which is obtained neglecting the quantum kinetic energy term in the equa- 
tion for the condensate, turns out to be very accurate |T^-|TP[. In this limit the Popov 
mean-field equations show an important scaling behaviour that we have recently investi- 
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gated in Ref. |3T|]. The scaling behaviour is best illustrated by introducing the dimen- 



sionless quantities pi = ^l/2mkBT^Pi, fi = ^m/2kBT^uJiri, jl = fi/kBT^, e = e/ksTl^, 
fiT = nT{2kBT^/mio'^Y^'^/N, = nQ{2kBT^/rmj'^Y^'^/N and the reduced temperature 
t = T/T^. By neglecting the quantum kinetic energy term in eq. ( ^2]) the mean-field 
equations can be rewritten in the following reduced form 

^o(^) = \(p--f'^- ^gnx) 9{fl - - 2ghT) , (68) 



r) = y^(p2 + r2 - + 2g{no + nT)f - g'^nl , (69) 



and 



/,p(_|)/(,7,) , (70) 



where the renormalized coupling constant g is given by 



^ = 87rr/5/Vl5~5.17(ArV6^) , (71) 



and is hence fixed by the interaction parameter rj defined in (^. In eq. (|7ID 6{x) is the step 
function and f{e/t) = (exp(e/t) — 1)^^ is the Bose function in terms of the reduced variables. 
The normalization condition for the reduced densities reads: / (if (no + nx) = 1. Equations 
(|68|)- ([70|) exhibit the anticipated scaling behaviour in the variables t and rj [^. Starting 
from their solutions one can calculate the condensate and the thermal densities as well as 
the excitation spectrum. This gives access to all the relevant thermodynamic quantities of 
the system. For example the condensate fraction is given by 



iV, 
N 

In particular the vanishing of the right hand side of eq. ( |72D fixes the value of the critical 
temperature tc = T^/T^. To the lowest order in the coupling constant g one finds |jl3| 



tc^l- 0.43r/^/2 = 1 - 1.3N^^^a/aHo ■ (73) 
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The energy of the system can be calculated starting from the standard thermodynamic 
relation in terms of the entropy (see eq. (0)). In units of NksTj^ one finds: 

E 5 /■* , , 5s 



7" + /o*'*'5^' 



where the entropy per particle s is given by (see eq. 

s(t,r/) = fcB^^^/ didp(^^f{~e/t)-\og{l-eM-m)) ■ (75) 

Finally, the ratio Q/Qrig of the moment of inertia to its rigid value (see eqs. (|5^), (|55|)) 
takes the form 

_e_ ^ 2 1 / dHpjx' + y')p'fm{l + fji/t) 

e„, 37r3c(3)t J df{x^ + r){Mf)+Mr)) ■ ^ ^ 

In the next sections we will show that the scaling behaviour is well achieved for the 
choice of parameters corresponding to the available experimental conditions. Scaling is 
consequently expected to provide a powerful tool for a systematic investigation of these 
trapped Bose gases. 

E. Elementary excitations and low temperature behaviour 

In this section we discuss the thermodynamic behaviour of a trapped Bose gas at low 
temperature (A;_bT <^ fi). Though at present this regime is not easily available in exper- 
iments, its theoretical investigation is highly interesting and reveals unexpected features 
concerning the competition between collective and single-particle effects which are worth 
discussing. 

We will derive our results in the framework of the scaling regime discussed in section 
2-D, where one assumes that the number of atoms in the trap is large enough to justify 
the use of the Thomas-Fermi approximation for the condensate and that the excited states 
can be described in the semiclassical approximation. This allows one to write the relevant 
thermodynamic functions in the form of integrals in coordinate and momentum space (see 
eqs. ([72|), ([74|)-([7^)) which are well suited for an analytic expansion at low temperature. 
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This expansion corresponds to exploring the region where ksT fi. At the same time 
the thermal energy ksT must be larger than the typical quanta of energy characterizing 
the elementary excitations of the system, in order to justify the use of the semiclassical 
approximation. In the absence of two-body interactions the quantum of energy is fixed by 
the oscillator energy hu. In the presence of interactions one has to distinguish between two 
regions. Inside the condensate the elementary excitations are of collective nature. For a 
spherical trap the dispersion law of these excitations can be expressed analytically in terms 
of the angular momentum quantum number i and of the number of radial nodes according 



to the law [^| 



e(n, i) = huj{2n'^ + 2nl + 3n + ^f/^ . (77) 

Their typical energy separation is still of the order of the oscillator energy huj. Near the 
classical boundary the relevant excitations are instead of single-particle type and in first 
approximation can be described by the Hartree-Fock Hamiltonian ( p6D which in the Thomas- 
Fermi limit, and for spherical traps takes the simplified form 

l^ + f-V-i^^l. (78) 



where R = \j'l^lmijj'^ is the classical boundary of the condensate. Actually the Hartree- 
Fock Hamiltonian (|78D gives correctly the energy of the elementary excitations only for 
r > R (see eq. (|35D). In fact for r < i? it ignores Bogoliubov effects. The potential 
term of the Hamiltonian ( [78D tends to confine the low energy single-particle states near the 
boundary. The typical energy gaps characterizing the low energy spectrum of (^) are of 
the order fiuj{R/ anoY^^ ~ huj{n/hujY^^ . In the large N limit this energy is larger than the 
oscillator energy, but still smaller than fi. As a consequence one can find a useful range of 
temperatures where to apply the low t expansion. In a homogeneous gas this region would 
correspond to the regime of collective excitations (phonons) which entirely determine the 
thermodynamics of the system. In a trapped Bose gas both collective and single-particle 
excitations taking place near the classical boundary can be important in the determination 
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of the low temperature behaviour. Actually, as it will be shown, for most thermodynamic 
quantities the low t behaviour is governed by single-particle excitations. 

Let us first consider the thermal energy of the system. At low temperature one can 
neglect the t-dependence in the interaction term characterizing the excitation spectrum (|69D 
and the integral appearing in (171) can be calculated explicitly giving the result 



E 5,1/- ,3_3. VF+WkP 



NkBTO 7' 7rX{3)Jf<^ exp(v^^T2p^/t) - 1 

TT'^ClS) Jf>V^ exp((p^ + - fioj/t) - 1 

where we have used the dimensionless variables of section 2-D with /io = yu(T = 0)/fcBT° = r) 
which at T = fixes the size of the condensate R = -/Jlq (see eq. (^8])). The first integral on 
the right hand side (r.h.s.) of eq. (|79|) gives the contribution to the energy from the region 
inside the condensate where the Bogoliubov effects take place, the second integral yields the 
contribution from the region outside the condensate where the spectrum is particle-like. By 



taking the linear phonon approximation e = p^flgriQ for the excitation spectrum one finds 
that the first integral on the r.h.s. of (^) has a divergent behaviour near the boundary 
i? as a consequence of the vanishing of the local sound velocity. This rules out the typical 
phonon-type behaviour for the temperature dependence of the energy. In fact a careful 
expansion of the integral ( |7DD at low t shows that the relevant region contributing to the 
integral is near the boundary. The expansion yields the law 

^ Ir^ + A^f/^ , (80) 



NkBT^ 7 

where A is a numerical factor. This behaviour points out the single-particle nature of the 
leading contribution to the thermal energy. In terms of the density of states one finds the 
result, valid for e ^ 0, (7(e) ~ e^^"^ which should be compared with the typical phonon law 
^(e) ~ e^. 

An opposite result is obtained if one investigates the low temperature behaviour of the 
thermal depletion (see eq. (|72D) 
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Nt If ,3~j3~ P +9^0 1 



J f<\/un 



If _ _ 1 

In fact in this case the particle distribution function F{p, r), whose integral gives the thermal 
depletion, is enhanced at small p by the phonon contribution (see eq. (|3)) which provides 
the leading effect in the integral (|HT|) according to the law: 



Eq. ( P^ provides the generalization of the Ferrel law holding for Bose superfiuids to 



a trapped Bose gas. In this case the single particle contribution arising from the boundary 
region is a higher order effect and behaves as t^^"^. 

In the scaling limit one can also obtain a useful estimate for the quantum depletion of 
the condensate. At T = the density of non-condensate atoms is given, in the semiclassical 
approximation, by the integral 6n{r) = J dp/{2TThYv^{p,r), where v^{p,r) is given in eq. 
(^). By employing the Thomas- Fermi approximation one finds the result 

Sn{r) = ^|§y(/^o - fY'Oih - f') , (83) 



already discussed in Ref. ||35|. It is worth noticing that the relative local quantum depletion 
is given by 5n{Y)/n{Y) = ^/?>{a?n{Y) /nY^'^ which coincides with the well known result for 
the quantum depletion of a homogeneous Bose gas @]. By integrating (|83D over space one 
gets the total number of quantum depleted atoms 

6 N r Tf^ 

— = / di5fi{f) = — ^ ~ 0.177^ . (84) 

N J ^ ^ 6y2C(3) 

The above estimate of the quantum depletion agrees reasonably well with the numerical 
results of Ref. ||T^ obtained from the full solution of the Bogoliubov equations. 

Let us finally discuss the low temperature behaviour of the moment of inertia. The 
normal component (^) gets a contribution from both the collective and single-particle ex- 
citations. One can easily show that the leading contribution to the integral ( [7BD at low 
temperature is given by single-particle excitations and behaves as 
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In this case the phonon contribution is a higher order effect and behaves as t^. 



III. RESULTS 

In this section we present the numerical results obtained by solving the equations of 
the Popov approximation both below and above the critical temperature. In section 3-A we 
discuss the temperature dependence of the condensate fraction and of the peak density in the 
center of the trap. Other structural properties such as the temperature dependence of the 
density profiles and of the mean square radii of the condensed and non-condensed part of the 
cloud are discussed. In section 3-B we present the results for the temperature dependence 
of the energy and of the specif heat Finally in section 3-C we discuss the temperature 
dependence of the moment of inertia. 

In the condensed phase, T < Tc, the self-consistent procedure is schematically divided 
in the following steps: 

i) the equation for the condensate wave-function 



2 



+ Kxt(r) - + g{no{r) + 2nr(r)) 

2m 



$(r) = (86) 



is solved, using the method described in Ref. [^], for the condensate density no(r) and the 
chemical potential /i, while keeping fixed the number of particles in the condensate, Nq{T), 
and the density nxir) of thermally excited atoms. 

ii) The condensate density and the chemical potential found above are used to calculate 
the excitation energies 

e(p, = I ( ^ + - + '^anir)) - g'nUv)] , (87) 



which yield a new non-condensate density tit through the equation (see eqs. (pSD, (p9|)) 
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iii) From the thermal density ht the total number of non-condensed atoms is obtained 
by direct integration 



and a new value for the number of atoms in the condensate is obtained from the normalization 
condition 



iv) With the new density ut and the new value for No{T) steps (i)-(iii) are repeated 
until convergence is reached. We stop the iterative procedure when successive iterations 
give values for the chemical potential and the condensate fraction, Nq{T)/N, which differ 
by less than a few parts in 10^. 

In the non-condensed phase, T > Tc, the condensate density uq^t) vanishes and the 
self-consistent procedure becomes easier because we do not need to solve the equation (pUf ) 
for the condensate. 

In the present work we have studied axially symmetric harmonic traps with the confining 
potential given by 



where Ur is the radial frequency in the x — y plane. Clearly the calculation can be also 
applied to different confining potentials Vext- Before performing the calculation described 
above one must specify the trap parameters (frequency uj^- and deformation A = uOz/oJr) as 
well as the parameters describing the system (mass m of the particles, scattering length a and 
total number of atoms A^). For a given temperature T the above self-consistent procedure 
gives directly the chemical potential, /x, the number of atoms in the condensate, Nq{T), the 
condensate and non-condensate densities no and ut- 

We have first checked the validity of the semiclassical approximation employed in the 
present work by comparing (see Fig. 2) our results for the condensate fraction with the 




(89) 



N = Nq{T) + Nt . 



(90) 



(91) 
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ones recently obtained in Ref. by solving the full Bogoliubov equations (p6D at finite 



temperature. The various parameter characterizing the trap and the gas are the same in 
the two calculations. Despite the relatively small value of (=2000) the agreement is good 
also at low T, revealing that no significant errors in the analysis are introduced by the use 
of the semiclassical approximation for the excited states. 

A. Condensate fraction and density distributions 



In Ref. ||T^ we have already studied the effect of the interatomic potential on the temper- 



ature dependence of the condensate fraction and on the critical temperature. We found that 
a repulsive interaction among the atoms enhances the thermal depletion of the condensate 
with respect to the prediction of the non-interacting model for the same number of trapped 
atoms. The effect results in a shift of the critical temperature towards lower temperatures. 
In Ref. [|13| we have also given an analytic expression for the shift 6Tc = T^ — of the crit- 
ical temperature from the prediction /csT° = huj(N/({3)y^^ of the non-interacting model 
holding in the large limit. The result, valid to lowest order in the scattering length and 
for large values of A^, reads 

^ ^ _0.73^^Ar-i/3 _ i.3^Ari/6 . (92) 

The first term on the right hand side (r.h.s.) of the above equation is independent of the 
interaction and gives the first correction to T° due to the finite number of atoms in the trap 



(see discussion in section 2-C and |^0|). The second term on the r.h.s. of eq. (|92D accounts 
for the effect of interaction; it can be either positive or negative, depending on the sign of a, 
increases with N and depends only on the geometric mean u through the harmonic oscillator 
length ajjo = (h/mcuY^'^. In the currently available experimental situations the shift due to 
interaction is the dominant effect as soon as the number of particles is of the order of 10^ 
for the JILA trap (a/a^o = 7.35 x IQ-^, X = VS) and of the order of 10^ for the MIT trap 
{a/ttHo = 2.55 X lO^'^, A = 18/320). The shift of the critical temperature due to interaction 
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has a quite simple physical interpretation: the presence of repulsive interactions among 
the atoms acts in reducing the density of particles in the center of the trap. Thus lower 
temperatures are needed to reach there the critical density for Bose-Einstein condensation. 
The opposite happens in the case of attractive interactions. 

In Fig. 3 we present results for the temperature dependence of the condensate fraction 
obtained with our self-consistent mean-field method, corresponding to = 5 x 10^ Rb atoms 
in the JILA trap, and to = 2.9 x 10^ Na atoms in the MIT trap. In the same figure we 
also plot the predictions for the same configurations obtained switching off the two-body 
interaction. The enhancement of the thermal depletion due to interaction is significant over 
the whole temperature range and the corresponding shift of the critical temperature is also 



clearly visible, being in quantitative agreement with the analytic estimate (p^). In Fig. 3 we 
also show the accuracy of the scaling behaviour discussed in section 2-D. Both configurations 
correspond to the value t] = 0.45 of the scaling parameter. The solid curve is the scaling 
function obtained from (^) with rj = 0.45, and the figure clearly shows that scaling is 
very well verified for these configurations and that hence extremely different experimental 
situations can give rise to the same thermodynamic behaviour. Only very close to the 
A^ = 5 X 10^ points exhibit small deviations from the scaling behaviour. In fact close to 
Tc the scaling law (|72D is approached more slowly with increasing A^. It is interesting to 
notice that finite size effects, responsible for the deviations from the scaling law 1 — in the 
non-interacting model (see eq. (Q)), are still visible for these configurations, whereas they 
are strongly quenched in the presence of the interaction. 

In Fig. 4 we present results for the condensate fraction Nq/N as a function of the reduced 
temperature t for three different values of the scaling parameter t], covering the presently 
available experimental conditions (see also Table I). In the figure the open diamonds with 
the error bars are the result of the Monte-Carlo simulation of Ref. which correspond 
to the value rj = 0.33 and which are in good agreement with our predictions. The dots are 
the experimental results of Ref. [0. In the experiments the number of particles A^ varies 
with t, with the value of t] ranging from 0.45 to 0.39. The experiments exhibit smaller 
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deviations from the non-interacting curve with respect to the ones predicted by theory. One 
should however keep in mind that the measured value of T is obtained from the velocity 
distribution of the thermal particles after expansion. The identification of this temperature 
with the one of the system before expansion ignores the interaction with the condensate 
which is expected to produce an acceleration of the thermal cloud. A preliminary estimate 
based on the calculation of the interaction energy between the thermally excited atoms and 
the condensate, shows that for T ~ 0.57° the final kinetic energy of the thermal cloud could 
be about 10% larger than its value before expansion. 

In Fig. 5 we show the temperature dependence of the density of particles in the center 
of the trap for two configurations: N — 10^ Na atoms in the MIT trap and N — 5000 Rb 
atoms in the JILA trap. The BEG transition is clearly marked by the discontinous change 
of slope in the two curves. 

In the limit of large numbers of particles the density in the center of the trap exhibits a 
scaling behaviour in terms of the parameters t and 77 when expressed in the dimensionless 
form discussed in section 2-D. In Fig. 6 the peak value of the reduced density n(f = 0) = 
n(r = 0){2kBT^/muj'^f/^/N is reported as a function of the reduced temperature t for three 
values of the scaling parameter rj (see also Table II). At T = one has h{r = 0) = Snrj'^^'^/lb 
which decreases by increasing rj. For T 3> T° the system approaches a classical ideal gas in 
the confining potential T4xt where h{r — 0) — (vri)"^/^. Notice that in the non- interacting 
model the peak density does not exhibit scaling behaviour for T < T°. For example at 
T = one has n(f = 0) = {2/nf/^{N/C{3)y/^ which diver ges as iV — >■ 00. 

In Figs. 7 a-d we show the density profiles of the total, condensate and non-condensate 
densities for four values of T/T° obtained from the self-consistent solution of the equations 
of the Popov approximation for the configuration corresponding to iV = 5000 Rb atoms 
in the JILA trap. The density profiles are plotted as a function of the radial coordinate. 
Similar profiles are obtained for the z-direction. At T/T° = 0.3 (7-a) the thermal density 
Ut is very small compared to the condensate density which almost coincides with the total 
density. By increasing T/T^ (7-b and 7-c) ut increases significantly and at T/T^ — 0.9 (7-c) 
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it becomes comparable with the condensate density in the center of the trap. Prom Figs. 
7-b, 7-c it clearly appears that the thermal density is not monotonically decreasing, and 
reaches a flat maximum near the boundary of the condensate. This behaviour arises from 
the repulsive interaction with the atoms in the condensate. Above T° (7-d) the condensate 
has disappeared and the thermal cloud has significantly spread out. In Figs. 7 a-d the scale 
on the vertical axis has been kept the same to give a visible evidence of the decrease of the 
density in the central region of the trap. The figure clearly shows the distinct behaviour 
of the two components of the gas: the condensate occupies the internal region of the trap 
while the thermal component forms a much broader cloud. 

In Fig. 8 we show the mean square radii of the condensate and of the thermal component 
as a function of T/T° for the same configuration of Fig. 7. The short-dashed hne refers to the 
mean square radius in the x-direction of the condensate atoms. By increasing T it decreases 
slowly and approaches a finite value at the transition point. The average radius of the 
thermal cloud is represented by the long-dashed line. It is always larger than the condensate 
radius and increases almost hnearly with temperature. The solid line corresponds to average 
radius of the total cloud / d^vn{v)x^ /N which coincides with the condensate radius at T = 
and goes over to the non-condensate radius for T > Tc. 

In the limit of large N the radii of the atomic cloud exhibit a scahng behaviour when 



expressed in the reduced units of section 2-D (fj = Jm/2kBT^uJiri). In Figs. 9-10 we 



show how the scaling behaviour is approached, along the x (y (x^)) and z-direction {y (z^)) 
respectively, for the JILA (open circles) and the MIT (solid circles) configuartions. Fig. 9 



shows that the scaling behaviour for ^J {x^) is well verified for both configurations, whereas 
along the z-direction the configuration with the smallest number of particles (N — 5 x 10^) 
exhibit a slight deviation from the scaling behaviour (see Fig. 10). This reveals that in the 
JILA trap finite size effects along the direction of stronger confinement have not completely 
disappeared for this value of A^. In the same figures we also plot the result for the non- 
interacting model [r] = 0). The dotted line is the result in the large A^ scaling limit, while 
the open and solid diamonds take into account finite size effects. The deviations from the 
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scaling law are much more evident in the non-interacting case. 

In Fig. 11 we plot the reduced mean radius of the cloud (ff) as a function of the reduced 
temperature t for three values of the scaling parameter rj (see also Table 111). The increase of 
the radius due to interaction is clearly visible by comparing the three curves with the result 
holding in the non-interacting model for t < 1: ^JJ^ = ^jE/GNksT^ = ^C(4)/2C(3)t2 
(see eq. (^)). Notice that for spherical traps the mean square radius is proportional to the 
harmonic potential energy {rf) = Euo/'iNkBT^- 

B. Energy and specific heat 

In this section we discuss the energetics of the trapped gases by presenting results for the 
chemical potential, total energy, release energy and specific heat. For all these quantities 
the scaling behaviour is well verified in the configurations realized in the experiments and 
we will therefore discuss our results only in the scaling limit (see section 2-D). 

In Fig. 12 we report the results for the chemical potential in units of fc^T^ as a function 
of the reduced temperature t for three values of the scaling parameter rj (see also Table 
IV). Notice that for t — the plotted quantity coincides with rj (see definition (|])). In the 
classical limit, T ^ T°, the dependence on the interaction parameter disappears and one 
finds the classical ideal gas prediction yu/fc^^c = ^log(C(3)/i'^)- 

In Fig. 13 we present the results for the total energy E of the system (see also Table V). 
At high temperature the behaviour is given by the classical law E/NksT^ = 3t. In the BEC 
phase the energy per particle is significantly larger than the prediction of the non-interacting 
model. This feature has been confirmed by the first experimental results on the temperature 
dependence of the release energy, which is the sum of the kinetic and interaction energy . 
In Fig. 13 the theoretical predictions for the release energy per particle, in units of ksT^, 
are shown together with the experimental results of Ref. (see also Table VI). Though 
the experimental points below Tc lie well above the non-interacting curve giving evidence 
for interaction effects, the accuracy is not enough to make a quantitative comparison with 
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theory. 

Finally in Fig. 15 we plot the specific heat per particle as a function of the reduced 
temperature t. These curves coincide with the derivative with respect to t of the curves of 
Fig. 13. The effect of interaction is to round off the peak at the transition with respect 
to the non-interacting model. Direct experimental results on the specific heat are not yet 
available. An indirect estimate can be obtained by differentiating the experimental curve of 
the release energy (see Ref. [H). One should however keep in mind that the release energy 
is roughly a factor two smaller than the total energy. 

C. Moment of inertia 

The superfluid properties of the trapped gases can be investigated by calculating the ratio 
between the moment of inertia of the system B and its rigid value Qrig (see section 2-A). 
In Fig. 16 our predictions for Q/Q^ig are shown as a function of the reduced temperature 
t. We find that, differently from the other quantities discussed above, this ratio does not 
exhibit a significant dependence on the scaling parameter t]. The figure also shows the very 
different behaviour exhibited by the non-interacting model where, even for large values of 
A^, there is a strong dependence on the number of trapped atoms and the anisotropy A of 
the trap (see eq. (|67|) ). 

IV. COMPARISON WITH THE HOMOGENEOUS BOSE GAS 

The thermodynamic behaviour of the trapped Bose gas discussed in the previous sections 
exhibits important differences with respect to the case of the uniform gas. The purpose of 
this section is to discuss these differences and to explain their physical origin. We will 
mainly focus on the temperature dependence of the condensate fraction where the effects 
are particularly visible and significant. 

Let us first discuss the behaviour of the ideal non-interacting gas. In the uniform case 
the critical temperature for the onset of Bose-Einstein condensation is given by the well 
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known expression 



2 / \ 2/3 



as a function of the density n of the gas. In the presence of harmonic confinement the 
critical temperature instead depends on the number of atoms and on the oscillator frequency 
according to the expression (H) 

ksT! = huj j , (94) 

holding in the large limit. By introducing the reduced temperature t = T/Tj? one can 
naturally compare the thermodynamic behaviour of the two systems. In particular the 
temperature dependence of the condensate fraction follows the laws 

^ = (95) 

and 

^ = 1-^', (96) 

in the uniform and in the harmonically trapped gases respectively. Different laws for the 
temperature dependence of the condensate are predicted to occur by changing the form of 
the trapping potential p. 

A major difference between the two systems concerns the relative shape of the conden- 
sate and of the thermally excited densities (see Fig. 7). In fact, differently from the uniform 
gas where the two components obvously overlap everywhere, in the presence of the external 
harmonic potential the size of the thermal component is larger than the one of the con- 
densate. The former is in fact of the order of the classical thermal radius {2kBT /muj'^Y^'^ ^ 
while the latter is fixed by the harmonic oscillator length ano = {h/rmuy^'^. As a con- 
sequence, as clearly shown by Fig. 7, most of the atoms of the thermal component lie 
outside the condensate and form a much more dilute gas, practically for any value of T. 
This behaviour has important consequences when we switch on the two-body interaction. 
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The effect of the repulsive interaction is in fact twofold. One the one hand it produces an 
expansion of the condensate whose size can significantly increase with respect to the one 
of the non- interacting model (this effect is fixed by the dimensionless parameter Na/ano 
as discussed in Refs. [p!"^-[T9|). On the other hand at finite temperature the repulsive forces 
favour the depletion of the condensate. In fact atoms are energetically favoured to leave the 
condensate where the repulsion is strong, and to reach the non-condensate phase which is 
much more dilute and interaction effects are consequently less important. This is exactly the 
opposite of what happens in the homogenous gas where the repulsion effects are stronger in 
the thermal component than in the condensate. This different behaviour can be understood 
by considering, for example, the Hartree-Fock expression for the interaction energy (last 
term of eq. (|58|), or, more directly, the Hartree-Fock expression for the particle distribution 
function (we use here the Hartree-Fock formalism because the effect can be discussed in a 
particularly transparent way in this approximation; analogous results are obtained in the 
Popov approximation as proven by the numerical results reported in Fig. 17 (a)), 

F(p, r) = (exp((pV2m + V,,t - + 2gn)/kBT) - l)"' , (97) 

which is obtained from eqs. ( p7[ ) and (^). In the homogeneous case one has Vext = 0, 
/i = 2gn — gnQ and the interaction effects produce a quenching of F and hence of the thermal 
depletion N-r = J dpdr/ {2TTh)^F{p, r) with respect to the ideal gas. In the case of harmonic 
confinement one sees that the relevant region of space contributing to Nj- lies outside the 
condensate where the term 2gn is negligible and interaction effects enter only through the 
chemical potential. In this case the repulsive interactions produce an enhancement of F with 
respect to the ideal gas, since the value of fi is positive for T < Tc. Notice that the value of 
the chemical potential is strongly affected by the interaction effects in the condensate and 
consequently the final effects on the thermal depletion can be sizable. In Fig. 17 (a-b) we 
show explicitly the temperature dependence of the condensate fraction in the uniform (a) 
and in the trapped (b) gases. In order to make the comparison significant, the density of 
the uniform gas has been fixed equal to the peak density n{r = 0) of the trapped gas at zero 
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temperature. The figure clearly shows that the effect of the interaction goes in the opposite 
direction in the two cases. 

An important consequence of the above results is that in a trapped Bose gas the repulsive 
interactions have the effect of lowering the critical temperature. This can be also understood 
by noting that the interaction produces a general expansion of the system and one has finally 
to deal with a more dilute gas with respect to the non-interacting gas with the same value 
of A^. In the uniform case, where the density is kept fixed, interaction effects are known 
to produce instead an increase of Tc, at least for very dilute gases. The corresponding 
shift cannot however be calculated in mean field theory, being associated with many-body 
fluctuation effects typical of critical phenomena. Also the behaviour of the condensate 
fraction near the critical point looks different in the two cases (see Fig. 17). In the uniform 
gas the mean field approach gives rise to a sizable gap in the condensate fraction at Tc. 
This gap, which is of the order of (a'^n)^/^, is an artifact of the theory and vanishes when 
fluctuation effects are taken into account near the critical point |^6[. This gap is present also 
in the case of the trapped gas, though it has a very small effect on the condensate fraction 
and on other thermodynamic functions. The reason is that the size of the condensate cloud 
near the transition point is small and the gap in the condensate fraction turns out to be of 



the order of r] |37 



Figure 17 points out also another important difference between the uniform and trapped 
gases. In the latter case phonons are always important in the calculation of the condensate 
fraction (this is proven by the significant difference between the prediction of the Hartree- 
Fock theory and the full Popov calculation which includes phonon effects). On the other 
hand in the case of harmonic confinement the Hartree-Fock approximation turns out to be 
an excellent approximation also at relatively low temperatures. In fact for these tempera- 
tures most of the contribution to thermodynamic properties of the trapped gas arise from 
excitations lying close to the classical boundary where the Hartree-Fock theory is a good 
approximation. 

It is finally worth noting that the ratio between the T = value of the chemical potential 
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and the critical temperature depends on the gas parameter a?n quite differently in the two 



the uniform gas the dependence is much stronger and is given by {a?nY/'^. 

In conclusion we have shown that a Bose gas in harmonic traps exhibits a thermodynamic 
behaviour quite different from the one of the uniform gas. The difference mainly originates 
from the different structure of the elementary excitations and from their coupling with the 
condensate. This should be taken into account when one compares the predictions of many- 
body theories obtained in uniform gases with the properties of these novel trapped systems. 



In this paper we have presented a systematic discussion of the thermodynamic behaviour 
of an interacting Bose gas confined in a harmonic trap. The formalism was based on an 
extension of Bogoliubov theory to non uniform systems at finite temperature. By using the 
semiclassical WKB approximation for the excited states we have been able to investigate 
various thermodynamic functions for a wide range of the relevant parameters of the system. 
The main results emerging from our analysis are briefly summarized below: 

i) For large numbers of atoms in the trap the effects of two-body interactions (limited in 
this work to the repulsive case) can be accounted for by the scaling parameter 



which gives the ratio between the T = value of the chemical potential and of the critical 
temperature for Bose-Einstein condensation. Typical values of r] in the available experiments 
range between 0.38 to 0.45. Physical systems with very different values of A^, scattering 
length and oscillator frequencies are predicted to give rise to the same thermodynamic 
behaviour provided they correspond to the same value of rj. 

ii) Interactions provide a significant quenching of the condensate fraction with respect 
to the non-interacting model and a shift of the critical temperature towards lower values. 



cases. In fact while in the trapped gas this ratio behaves as {a?n{v = 0)Y^^ (see eq. (|)), in 



V. CONCLUSIONS 




(98) 
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This behaviour, which differs from the one exhibited by a uniform Bose gas, points out 
a unique and interesting feature exhibited by these trapped systems: the occurrence of 
elementary excitations located outside the condensate. The thermal excitation of these 
modes is energetically favoured by the presence of the repulsive interaction and plays a 
crucial role in the thermodynamics of the system. 

iii) Interactions are also responsible for a significant enhancement of the thermal exci- 
tation energy. Evidence for such an effect has been recently become available from first 
measurements of the release energy. 

iv) Results have been also obtained for the moment of inertia whose deviations from the 
rigid value point out the occurrence of superfluidity. Such deviations are significant in a 
rather wide range of temperatures below and have been found to depend very little on 
the value of the interaction parameter rj. 

v) Explicit comparison with the predictions of Path Integral Monte Carlo simulations 
as well as with the full solutions of the discretized Bogoliubov equations (without using the 
WKB approximation) have revealed that the method presented in this paper is remarkably 
accurate in the relevant range of temperatures. Experimental results instead are not yet 
enough accurate to check the theoretical predictions in a quantitative way. 



ACKNOLEDGMENTS 



We gratefully acknowledge useful discussions with Franco Dalfovo and Allan Griffin. We 



would also like to thank E. Zaremba for providing the numerical results of Ref. |14 



38 



REFERENCES 



[1] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman and E.A. Cornell, Science 
269, 198 (1995); C.C. Bradley, C.A. Sackett, J.J. Toilet and R.G. Hulet, Phys. Rev. 
Lett. 75, 1687 (1995); K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van Druten, D.S. 
Durfee, D.M. Kurn and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 

[2] M.-O. Mewes, M.R. Andrews, N.J. van Druten, D.M. Kurn, D.S. Durfee and W. Ket- 
terle, Phys. Rev. Lett. 77, 416 (1996). 

[3] J.R. Ensher, D.S. Jin, M.R. Matthews, C.E. Wieman and E.A. Cornell, Phys. Rev. 
Lett. 77, 4984 (1996). 

[4] N. Bogoliubov, J. Phys. USSR 11, 23 (1947). 

[5] K. Huang and C.N. Yang, Phys. Rev. 105, 767 (1957); K. Huang, C.N. Yang and J.M. 
Luttinger, Phys. Rev. 105, 776 (1957); T.D. Lee and C.N. Yang, Phys. Rev. 105, 1119 
(1957); T.D. Lee, K. Huang and C.N. Yang, Phys. Rev. 106, 1135 (1957). 

[6] E.P. Gross, Nuovo Cimento 20, 454 (1961); L.P. Pitaevskii, Zh. Eksp. Teor. Fiz. 40, 
646 (1961) [Sov. Phys. JETP 13, 451 (1961)]. 

[7] A.L. Fetter and J.D. Walecka, Quantum Theory of Many Particle Systems (McGraw- 
Hill, New York, 1971); A.L. Fetter, Ann. Phys. (N.Y.) 70, 67 (1972); A.L. Fetter, Phys. 
Rev. A 53, 4245 (1996). 

[8] S.R. de Groot, G.J. Hooyman and C.A. Ten Seldam, Proc. R. Soc. (London), A 203, 
266 (1950); V. Bagnato, D.E. Pritchard and D. Kleppner, Phys. Rev. A 35, 4354 (1987). 

[9] V.V. Goldman, I.F. Silvera and A.J. Leggett, Phys. Rev. B 24, 2870 (1981); D.A. Huse 
and E.D. Siggia, J. Low Temp. Phys. 46, 137 (1982). 

[10] J. Oliva, Phys. Rev. B 39, 4197 (1989). 

[11] T.T. Chou, C.N. Yang and L.H. Yu, Phys. Rev. A 53, 4257 (1996); T.T. Chou, C.N. 

39 



Yang and L.H. Yu, |cond-mat /9605058 . 



[12] A. Griffin, Phys. Rev. B 53, 9341 (1996). 

[13] S. Giorgini, L.P. Pitaevskii and S. Stringari, Phys. Rev. A 54, 4633 (1996). 
[14] D.A.W. Hutchinson, E. Zaremba and A. Griffin, |cond-mat/ 96 11023, 



[15] A. Minguzzi, S. Conti and M.P. Tosi, Jou. of Phys. Cond. Matter, 9, L33 (1997). 
[16] W. Krauth, Phys. Rev. Lett. 77, 3695 (1996). 

[17] M. Edwards and K. Burnett, Phys. Rev. A 51, 1382 (1995); P.A. Ruprecht, M.J. Hol- 
land, K. Burnett and M. Edwards, Phys. Rev A 51, 4704 (1995). 

[18] G. Baym and C. Pethick, Phys. Rev. Lett. 76, 6 (1996). 

[19] F. Dalfovo and S. Stringari, Phys. Rev. A 53, 2477 (1996). 

[20] Corrections to (H) due to finite size effects have been the object of several recent papers 



[21] Actually the applicability of the semiclassical approximation in the region near the 
boundary of the condensate implies a more severe condition (see Sect. 2-E). 

[22] D.S. Jin, J.R. Ensher, M.R. Matthews, C.E. Wieman and E.A. Cornell, Phys. Rev. 
Lett. 77, 420 (1996); D.S. Jin, M.R. Matthews, J.R. Ensher, C.E. Wieman and E.A. 
Cornell, preprint. 

[23] M.-O.Mewes, M.R.Andrews, N.J. van Druten, D.M. Kurn, D.S. Durfee and W. Ketterle, 
Phys. Rev. Lett. 77, 988 (1996). 

[24] K.G. Singh and D.S.Rokhsar, Phys. Rev. Lett. 77, 1667 (1996); S. Stringari, Phys. Rev. 
Lett. 77, 2360 (1996). 

[25] M. Edwards, P.A. Ruprecht, K. Burnett, R.J. Dodd and C.W. Clark, Phys. Rev. Lett. 
77, 1671 (1996). 

40 



[26] P.C. Hohenberg and P.C. Martin, Ann. Phys. (N.Y.) 34, 291 (1965). 

[27] V.N. Popov, Functional integrals and collective excitations (Cambridge University Press, 
Cambridge, 1987). 

[28] See for example L.D. Landau and E.M. Lifshitz Statistical Physics (Pergamon, Oxford, 
1980) Part 1. 

[29] S. Stringari, Phys. Rev. Lett. 76, 1405 (1996). 

[30] S. Grossmann and M. Holthaus, Phys. Lett. A208, 188 (1995) and Zeit. f. Naturforsch. 
50a, 323 (1995); W. Ketterle and N.J. van Druten, Phys. Rev. A 54, 656 (1996); K. 
Kirsten and D.J. Toms, Phys. Rev. A 54, 4188 (1996); H. Haugerud, T. Haugset and 
F. Ravndal, Phys. Lett. A 225, 18 (1997). 

[31] S. Giorgini, L.P. Pitaevskii and S. Stringari, submitted to Phys. Rev. Lett. . 

[32] In the hmiting case 1] ^ the equations can be solved analytically. For example for the 
condensate fraction (see eq. (173)) one finds the expansion: Nq/N = {l—t^)—({2)rit'^{l — 



[33] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

[34] R.A. Ferrel, N. Meynhard, H. Schmidt, F. Schwabl and P. Szepfalusy, Ann. Phys. (N.Y.) 
47, 565 (1968). 



[35] E. Timmermans, P. Tommasini and K. Huang, |cond-mat/ 9609234 . 



[36] M Bijlsma and H.T.C. Stoof, Phys. Rev. A 54, 5085 (1996). 

[37] Notice that the same type of difficulties arise also for the space distribution no(r) of the 
condensate density in the Thomas- Fermi limit (see eq. (§^). The reason is that near the 
boundary critical fluctuations are important not only near the transition temperature. 
See K. Damle, T. Senthil, S.N. Majumdar and S. Sachdev, Europhys. Lett. 36, 7 (1996). 



41 



FIGURES 

FIG. 1. Temperature dependence of the condensate fraction in the non-interacting limit. The 
sohd circles correspond to the exact quantum calculation for N = 1000 atoms in the JILA trap 
and the solid line to the semiclassical approximation of eq. (|64|). The short-dashed line refers to 
the large limit. 

FIG. 2. Condensate fraction as a function of T/T^. The solid circles are the result of the calcula- 
tion of Ref. [^] for N = 2000 Rb atoms in an isotropic trap with oscillator length auo = 7.62 x 10~^ 
cm, obtained without using the semiclassical approximation. The solid line is the result for the 
same configuration employing the semiclassical approximation. Also plotted is the curve of the 
non- interacting model in the large N limit (short-dashed line). 

FIG. 3. Theoretical predictions for the condensate fraction as a function of T/T^. The open 
circles correspond to A^ = 5 x 10^ Rb atoms in the JILA trap {a/ auo = 7.35 x 10^'^, A = \/8). 
The solid circles correspond to A^ = 2.9 x 10'' Na atoms in the MIT trap {a/auo = 2.55 x 10^'^, 
A = 18/320). The solid line corresponds to the scaling limit with r/ = 0.45. The dotted line is 
the 1 — curve of the non-interacting model in the large A^ limit. The open and solid diamonds 
correspond to A^ = 5 x 10^ and A^ = 2.9 x 10^ non-interacting particles in the JILA and MIT traps 
respectively. 

FIG. 4. Condensate fraction in the scaling limit as a function of T/T^ for three values of rj. 
Solid line: rj = 0.45, long-dashed line: rj = 0.39, short-dashed line: rj = 0.31. Open diamonds: 
PIMC results of Ref. ||l^. Solid circles: experimental results from Ref. [^. The dotted line refers 
to the non- interacting model in the large A^ limit. 

FIG. 5. Peak density as a function of T/T^. Solid line: A^ = 10^ Na atoms in the MIT trap, 
long-dashed line: N = 5 x 10^ Rb atoms in the JILA trap. 

FIG. 6. Reduced peak density as a function of T/T^ for three values of the scaling parameter 
r/ (see Fig. 4). 
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FIG. 7. Density profiles as a function of the radial coordinate at z = for N = 5x 10^ Rb atoms 
in the JILA trap at different temperatures. Solid line: total density, long-dashed line: condensate 
density, dotted line: non-condensate density. Lengths are in units of the radial oscillator length 
Qr = {fi/mujr)^/'^ . Densities are in units of a~^. 

FIG. 8. Root mean square (r.m.s.) radii, in units of a^, in the x-direction for N = 5 x 10^ 
Rb atoms in the JILA trap as a function of T/T^. Solid line: average radius of the total cloud; 
long-dashed line: average radius of the thermal cloud; short-dashed line: average radius of the 
condensate cloud. 

FIG. 9. r.m.s. radii in the x-direction in reduced units as a function ofT/T^ for rj = 0.45 (solid 
line). The open circles refer to = 5 x 10^ Rb atoms in the JILA trap. The solid circles correspond 
to A = 2.9 X 10'' Na atoms in the MIT trap. The dotted line refers to the non-interacting model 
in the large A limit. The open and solid diamonds correspond to A = 5 x 10^ and A = 2.9 x 10^ 
non-interacting particles in the JILA and MIT traps respectively. 

FIG. 10. r.m.s. radii in the z-direction (see Fig. 9). 

FIG. 11. r.m.s. radii in reduced units as a function of T/T^ for three values of the scaling 
parameter rj (see Fig. 4). The dotted line refers to the non-interacting model in the large A limit. 

FIG. 12. Chemical potential as a function of T/T^ for three values of the scaling parameter rj 
(see Fig. 4). The dotted line refers to the non-interacting model in the large A limit. 

FIG. 13. Total energy of the system as a function of T/T^ for three values of the scaling 
parameter rj (see Fig. 4). The dotted line refers to the non-interacting model in the large A limit. 

FIG. 14. Release energy of the system as a function of T/T^ for three values of the scaling 
parameter rj (see Fig. 4) . The solid circles are the experimental results of Ref . . The dotted line 
refers to the non-interacting model in the large A limit. 
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FIG. 15. Specific heat as a function of T/T^ for three values of the scahng parameter r] (see 
Fig. 4). The dotted fine refers to the non- interacting model in the large N limit. 

FIG. 16. Ratio Q/Qrigid as a function of T/T'^ for three values of the scaling parameter 77 (see 
Fig. 4) . The three curves coincide almost exactly and are represented by the solid line. The dotted 
line refers to iV = 5 x 10^ atoms in the JILA trap in the non-interacting model. The dot-dashed 
line refers to N = 2.9 x 10^ atoms in the MIT trap in the non-interacting model. 

FIG. 17. Condensate fraction as a function of T/T°. (a) refers to a homogeneous gas with 
the gas parameter a^n = 6.2 x 10~^, (b) refers to a trapped gas with the same value of the gas 
parameter in the center of the trap at T = 0, corresponding to the value rj = 0.45. Solid line: 
Popov approximation; long-dashed line: Hartree-Fock approximation; dotted line: 1 — t^^^ (a) and 
1 — (b) curves of the non-interacting model in the homogeneous and trapped case respectively. 
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TABLES 



TABLE I. Condensate fraction 





No/N 


No/N 


No/N 




r] = 0.31 


T] = 0.39 


r] = 0.45 


0.09 


0.99 


0.99 


0.99 


0.25 


0.95 


0.94 


0.94 


0.50 


0.76 


0.73 


0.71 


0.75 


0.38 


0.33 


0.29 


0.91 


0.08 


0.05 


0.03 




TABLE II. Peak density in 


reduced units 






n(f = 0) 


h{r = 0) 


n(f = 0) 




rj = 0.31 


7/ = 0.39 


7/ = 0.45 


0.25 


3.40 


2.43 


1.95 


0.50 


3.19 


2.28 


1.81 


0.75 


2.55 


1.79 


1.41 


1.00 


0.30 


0.28 


0.26 


1.25 


0.15 


0.15 


0.15 
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TABLE III. Root mean square radius in reduced units 





/ / ~2\ 


/ / ~2\ 


/ / ~2\ 




rj = 0.31 


rj = 0.39 


Tj = 0.45 


0.25 


0.22 


0.24 


0.26 


0.50 


0.29 


0.31 


0.33 


0.75 


0.46 


0.48 


0.50 


1.00 


0.67 


0.68 


0.68 


1.25 


0.77 


0.78 


0.78 



TABLE IV. Chemical i)<)(('iiiial 













r] = 0.31 


rj = 0.39 


T] = 0.45 


0.25 


0.31 


0.39 


0.45 


0.50 


0.30 


0.37 


0.43 


0.75 


0.24 


0.31 


0.36 


1.00 


0.02 


0.03 


0.04 


1.25 


-0.70 


-0.69 


-0.69 



TABLE V. Total energy 



T/ro 


E/NHbT^ 




E/NkBT^ 




ri = 0.31 


rj = 0.39 


r} = 0.45 


0.25 


0.24 


0.30 


0.35 


0.50 


0.47 


0.55 


0.61 


0.75 


1.27 


1.38 


1.47 


1.00 


2.73 


2.75 


2.77 


1.25 


5.59 


3.61 


3.62 
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TABLE VI. Release energy 





ER/NksT^ 


T — 1 / TV T 7 mC\ 


7^ / TV T 7 /~nC\ 

ER/NksT^ 




T] = 0.31 


J] = 0.39 


T] = 0.45 


0.25 


0.10 


0.12 


0.14 


0.50 


0.22 


0.25 


0.28 


0.75 


0.63 


0.68 


0.72 


1.00 


1.36 


1.37 


1.38 


1.25 


1.79 


1.80 


1.81 
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